r/o  u/i 


kO-AO»S  010  01SC0NSIN  UMVHUDIMN  MTMDWTICS  DC 
THK  MUmOM  TO  *  SMOOTH  TOC  CAM  K 

nm  oo  sc  ouch 

UNCLASSIFIED  Wtc-TM  SSSS 


in  ciii.iu) 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


ADA083818 


Mathematics  Research  Center  u 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 

February  1980 


/ 


(Received  November  7,  1979) 


Ap  prevail  far  paklic  release 
Distrikatiaa  eelimitatf 


iXrf  U.S.  Research  Office 

P.O.  Box  12211 
cz  Research  Triangle  Park 
North  Carolina  27709 


National  Science  Foundation 
Washington,  D.  C.  20550 


80  4  9  047 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 

THE  SOLUTIONS  TO  A  SMOOTH  PDE  CAN  BE  DENSE  IN  C[I] 

R.  C.  Buck 

Technical  Summary  Report  #2036 
February  1980 

ABSTRACT 

The  solutions  to  a  partial  differential  equation  that  arises  in  the 
physical  sciences  are  expected  to  be  restricted  in  nature,  since  they  are 
intended  to  represent  some  physically  significant  behavior.  In  particular, 
they  can  approximate  the  general  continuous  function  on  a  compact  set  K 
only  if  K  is  "thin",  e.g.  nowhere  dense  in  Rn.  The  purpose  of  this 
note  is  to  show  that  this  does  not  hold  for  all  smooth  partial  differential 
equations.  Specifically,  for  any  n  >  2  there  exist  partial  differential 
equations  of  polynomial  type  on  Rn  whose  C°°  solutions  are  uniformly 
dense  in  the  space  C[I]  of  all  continuous  functions  on  the  n-cube  I. 
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SIGNIFICANCE  AND  EXPLANATION 


The  solutions  to  a  partial  differential  equation  that  arises  in  the 
physical  sciences  are  expected  to  be  restricted  in  nature,  since  they  are 
intended  to  represent  some  physically  significant  behavior.  In  particular, 
they  can  approximate  the  general  continuous  function  on  a  compact  set  K 
only  if  K  is  "thin",  e.g.  nowhere  dense  in  Rn.  The  purpose  of  this  note 
is  to  show  that  this  does  not  hold  for  all  smooth  partial  differential 
equations.  Specifically,  for  any  n  >  2  there  exist  partial  differential 
equations  of  polynomial  type  on  R  whose  C  solutions  are  uniformly 
dense  in  the  space  C[I]  of  all  continuous  functions  on  the  n-cube  I. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


THE  SOLUTIONS  TO  A  SMOOTH  PDE  CAN  BE  DENSE  IN  C II] 
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The  solutions  to  a  partial  differential  equation  that  arises  in  the  physical 
sciences  are  expected  to  be  restricted  in  nature,  since  they  are  intended  to 
represent  some  physically  significant  behavior.  In  particular,  they  can  approximate 
the  general  continuous  function  on  a  compact  set  K  only  if  K  is  "thin”,  e.g. 
nowhere  dense  in  Rn.  The  purpose  of  this  note  is  to  show  that  this  does  not  hold 
for  all  smooth  partial  differential  equations.  Specifically,  for  any  n  >  2  there 
exist  partial  differential  equations  of  polynomial  type  on  R°  whose  C  solutions 
are  uniformly  dense  in  the  space  C[I]  of  all  continuous  functions  on  the  n-cube 


I. 


Here,  I  is  the  set  of  all  x  = 
The  vector  exponent  a  =  (a^,  a2#  ... 
j '■  |  =  J  ;  the  general  multinomial, 


(x, ,  x,,  ...  ,  x  )  E  Rn  with 
12  n 

,  a  >  has  integer  components 
n 

denoted  by 


0  <  Xj,  <  1. 

>  0,  and 


is  said  to  have  (total)  degree  |a|. 

We  also  use  this  notation  for  partial  derivatives,  so  that 

,  a  n  a. 

(t)  F=  n  {■£■)  ‘r 


i  3xi 


is  one  of  the  derivatives  of  F  of  order  |a|. 
We  prove  the  following  result: 


Main  Theorem.  If  n  >  2  there  are  integers  k  =  k(n)  and  A  »  A(n)  and  a  partial 
differential  equation  on  Rn  of  order  not  greater  than  k,  whose  C  solu¬ 
tions^  F(x)  are  uniformly  dense  in  C[I),  the  space  of  all  continuous  real 
valued  functions  on  the  n-cube  I.  The  form  of  the  differential  equation  is 

<1>  Q<V  w2,  ...  ,  wrj)  =0 
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<  k,  and  where  Q  is  a  nontrivial  polynomial  of  total  degree  at  most 


The  construction  of  such  a  "universal*'  PDE  which  we  will  outline  depends  upon 
a  combination  of  familiar  facts,  of  which  the  central  one  is  the  Kolmogorov  solu¬ 
tion  of  Hilbert's  13th  Problem  (See  (3]).  This  asserts  that  every  continuous 
function  F  on  I  can  be  represented  there  in  the  form 

2n+l 

(2)  F  =  l  f .  o  g. 


where  the  f ^  are  continuous  functions  of  one  variable,  and  the  g^  are  sums  of 
such  functions: 


(3) 


g^(x)  =  j  (xx  > 


.  i,  (x  ) 
3  2  2 


. :  ( x  ) 
:  n  n 


The  various  proofs  of  this  surprising  fact  depend  on  the  Baire  category 
theorem,  and  are  not  constructive.  It  is  known  however  that  one  cannot  if.  general 
hope  for  all  the  component  functions  f^  and  to  be  differentiable,  even  if 

F  itself  is  reasonably  smooth  (See  again  [3]). 

We  introduce  the  special  term  smooth  Kolmogorov  function  for  a  fun  -t ion  f 
defined  on  I  that  is  represented  there  in  the  format  (2)  and  (3),  where  each  of 
the  component  functions  is  of  class  C  .  As  observed  above,  these  for",  a  proper 
subset  of  CIX]  . 


Lemma  1 .  The  clans  of  smooth  Kolmogorov  functions  are  uniformly  dense  m  Ctrl. 

This  depends  on  the  elementary  observation  that  if  9n  !  A  •  B  and 

f  :  B  C,  are  continuous  mappings  between  compact  metric  spaces,  and  if  •  f n 

and  {g  }  converge  uniformly  on  their  domains  to  f  and  g,  respectively,  then 

f  o  g  converges  to  fog,  uniformly  on  A. 
n  n 

The  proof  of  the  main  theorem  will  be  completed  by  showing  that  a  PDE,  of  the 


form  (1) ,  can  be  constructed  whose  solutions  include  all  the  smooth  Kolmogorov 
functions.  The  argument  uses  two  elementary  results  dealing  with  polynomials. 


Lemma  2 .  The  collection  of  polynomials  in  r  indeterminates  of  total  degree  at 

fr  +  d-j 

most  d  forms  a  vector  space  V ( r , d )  of  dimension 


1  r  ) 


The  second  result  is  the  standard  theorem  on  algebraic  dependence. 


Lemma  3.  (i)  If  N  >m  and  P, (u) ,  P_(u),  ...  ,  p  (u)  is  a  collection  of  poly- 

■  -  -  — — -  12  N 

nomials  in  the  indeterminates  u  =  (u. ,  u_  ,  ...  ,  u  ),  then  there  is  a  non- 
- —  12  m 

trivial  polynomial  Q  such  that 

(4)  Q(P,,  P.,,  ...  ,  P>  =  0  . 

1  J  N 

(ii)  If_  deg  (P  )  _  for  each  i,  then 


(5) 


deg  (Q)  l  N 


m/N-m) 


Proof : 


Set  w.  =  P.(u,,  u„,  ...  ,  u  ),  for  i  =  1,  2,  ...  ,  N.  The  goal  is  to 
i  i  1  2  m 

eliminate  the  indeterminates  u^,  obtaining  (4).  Choosing  an  as  yet  unspecified 


form  all  the  multinomials  w  for 


By  Lemma  2,  we  will  have 


(6) 


s('.l  = 


of  these.  But, 
(7) 


w  =  P  <u  ,  u  ,  ...  ,  u  ) 
«ii2  m 

i=l 


which  is  a  polynomial  in  the  u^  of  total  degree  at  most  )  "  v  =  ■  ]  i|  _ 
belonging  to  the  vector  space  V(m,  ■>/,).  By  Lemma  2,  the  dimension  of  this  space  is 

"L\  +  m-i 


(0) 


fL,  ♦  m, 

t  (A)  = 

l  m  ) 


Thus,  if  s()  t(),  the  polynomials  will  be  linearly  dependent,  and  there  will 

»  0.J 

exist  scalars  c  ,  not  all  zero,  such  that  o  =  )  c  w‘=Tc  ::  w. 

«  i  L  a  ,  i 

i  a  1 

Q(w, ,  ...  ,  w  ),  a  polynomial  of  degree  at  most  such  that  (4)  holds.  To  obtain 
1 


-3' 


N 

s(A)  >  t(A),  we  observe  that  (6)  shows  that  s(A)  =  0(A  ),  while  (0)  gives 
t  ( A)  =  0(Am);  since  N  >  m,  s(A)  >  t(A)  when  A  is  sufficiently  large. 
The  explicit  estimate  given  in  (5)  is  easily  obtained,  since 


N-m 


The  proof  of  the  main  theorem  is  now  completed  by  returning  to  (2)  and  dif¬ 
ferentiating  it  repeatedly,  up  to  partial  derivatives  of  total  order  at  nost  k, 
thereby  obtaining  a  large  number  of  equations  of  the  form 


(9) 


(  1  F  =  F  *  P  (u  ,  u  ,  ...  ,  u 

dx  J  a  a  i-  2  r 


where  the  u^  are  the  functions  arising  from  differentiating  the  component  func¬ 
tions  entering  into  the  Kolmogorov  ^ori-iat  representation  of  F,  and  p  is  a 

a 

polynomial. 

Consideration  of  a  related  case  may  clarify  this  (See  11]).  Suppose  that 
(10)  G  (x^.Xj )  =  f  ( <x2 >  >  =  (f  °  9)  (x) 


where  f,  ©  and  v  are  in  C  .  By  differentiation,  we  obtain 


1.0 


U1  U2 


0,1 


u  u 
1  3 


wliere 


2,0 


1,1 


fog 
f”  o  g 


U,  uc  +  u, (u_ ) 
15  4  2 


u  u  u 
4  2  3 


etc 


tt'" 


etc. 


-4- 


In  this  example,  it  may  be  seen  that  the  can  be  eliminated,  resulting  in  a  poly¬ 

nomial  relation  among  the  G^,  namely: 

<G10)2  G01  G12  ‘  G10(G01)2  G21  '  <G10)2  G11  G02  +  (G01)2  G11  G20  =  ° 


This  is  a  partial  differential  equation  of  order  3  and  degree  4  satisfied  by  all  the 
functions  G  that  can  be  represented  in  the  format  (10)  with  smooth  component 
functions  f,  <t> ,  and  tji.  Note  that  this  PDE  does  not  depend  upon  f,  i,  or  . , 
but  only  upon  the  specific  format  displayed  in  (10) . 

The  procedure  in  the  Kolmogorov  case  is  essentially  the  same.  By  differentiating 
(2)  up  to  order  k,  we  produce  N  equations  of  the  form 

1  <  |a|  <  k 


r  =  p  (u  ,  u  ,  . . 
a  a  1  2 


V 

m 


where  the  u.  are  the  functions  rP  .  0.  and  (D^  f.)  o  g.  for  p 
i  It  3  1 


;  1,  2,  ...  ,  k, 

j  =  1,  2,  ...  ,  2n  +  1,  i  =  1,  2,  ...  ,  n,  and  D  is  d/dt.  Thus, 
m  »  k(2n  +  1) (n  +  1).  By  Lemma  2,  the  number  of  equations  obtained,  N,  is  the 
dimension  of  V(n,k)  minus  1,  since  we  omit  the  0th  order  derivative.  Hence, 

k  +  ni 


N 


-  1 


Moreover,  each  of  the  polynomials  has  degree  at  most  k  +  1,  and  integral 

coefficients.  Since  N  =  0(kn)  while  n  >  2  and  m  =  0(k),  we  see  that  N  >  m 
when  k  is  sufficiently  large.  Thus,  by  Lemma  3,  there  is  a  polynomial  Q  with 
rational  coefficients  and  of  degree  less  than  the  integer  A  given  in  (5)  such 
that 

Q<  •••  .  (  3^  )  F,  ...  )  =  0 

which  depends  only  upon  the  integer  n  and  the  format  given  in  (2)  and  (3)  and  not 
upon  the  specific  functions  f^  and  ^0^  that  appear  there.  Thus,  this  partial 
differential  equation  will  have  as  solutions  all  the  smooth  Kolmogorov  functions  of 


n  variables  on  I.  The  main  theorem  follows. 


Explicit  values  of  k  and  &  can  be  found  for  any  n  >2.  However,  it  is  no 
likely  that  anyone  will  wish  to  exhibit  Q  itself,  since  the  order  and  degree  for 
low  dimensions  n  are  so  large.  Interest  in  this  theorem  lies  in  the  fact  that 
such  "universal"  differential  equations  exist,  and  that  they  are  themselves  smooth 
(of  polynomial  type)  and  that  it  is  the  C  solutions  of  these  equations  that  have 
the  approximation  property. 

For  the  record,  we  record  that  if  n  =  2,  k  =  28,  tj  =  434,  m  =  420;  if  n  =  1 
k  =  4,  N  =  1000,  m  =  924;  that  if  n  =  32,  k  =  3,  N  =  6544,  m  =  6435.  All  of  the 
upper  bound  estimates  for  £  are  very  large. 

There  is  no  reason  to  believe  that  "universal"  differential  equations  of  much 
lower  order,  degree,  and  dimension  may  not  exist.  However,  the  present  line  of 
argument  does  not  yield  any  PDE  of  order  2.  It  would  be  of  interest  to  show  that 
no  second  order  partial  differential  equation  of  polynomial  type  in  n  variables 
can  have  the  approximation  property.  It  is  very  plausible  that  this  is  true  for 
polynomial  equations  in  one  variable,  of  any  order. 

The  restriction  to  equations  of  polynomial  type  may  be  an  essential  feature  of 

the  example  given  in  this  paper.  In  subsequent  discussions,  Professors  Crandall 

and  Turner,  (2] ,  have  shown  that  it  is  relatively  easy  to  construct  DDEs  of  the 

form  Q(u',u")  =  0  with  solution  sets  that  are  uniformly  dense  on  intervals,  and 
■  00 

with  Q  of  class  C  .  These  are  constructed  by  careful  selection  of  the  desired 
solution  sets  S  in  advance ,  and  then  construction  of  the  function  2  so  that  it 
vanishes  on  a  countable  family  of  curves  where 

y  (t)  =  (u'  (t)  ,u"  (t) ) ,  for  t  c  I  and  u  e  S  c  S  . 
u  0 

This  of  course  forces  SQ  to  be  a  subset  of  the  solutions  of  Q(u',u")  =  0.  In 

particular,  such  a  construction  of  SQ  is  possible  for  any  c”  function  Q  with 


compact  support. 
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